typeset using JPSJ.sty <ver.l.Ob> 



Spin Splitting in de Haas-van Alphen Oscillation in Two-Dimensional 

Two-Band Systems 

Keita Kishigi, Yasumasa Hasegawa and Mitake Miyazaki 

Faculty of Science, Himeji Institute of Technology, Akou-gun, Hyogo 678-1297, Japan 
(Received February 1, 2008) 

We study the effects of the Zeeman term on the de Haas van Alphen (dHvA) oscillation in 
two-dimensional two-band systems. We found that the Fourier transform amplitudes of the 
oscillations are not described by the spin reduction factor in the Lifshitz-Kosevich formalism in 
two-dimensional systems. The anomalous dependence on the effective p-factor can be observed 
by tilting-angle dependence of the dHvA oscillation in quasi-two-dimensional organic conductors 
and Sr2RuC>4. 
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The magnetization and magnetoresistance oscillate as 
a function of the inverse of the magnetic field (H) with 
the period proportional to the extreme area (/) of the 
closed orbit on the Fermi surface. These phenomena are 
known as the Shubnikov-de Haas (SdH) oscillation and 
the de Haas-van Alphen (dHvA) oscillation. The dHvA 
oscillation can be described in the seraidassical theory 
as Lifshitz and Kosevich (LK) formulaBB) 

In this study, we consider the system at the clean limit 
and zero temperature. In the two-dimensional single- 
band system, the magnetization, ^f^(n, H), for the fixed 
chemical potential (/z) is given by 
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where R£ is the reduction factor due to electron spin, 
g ~ 2 is the g-factor, to is the cyclotron effective mass 
and mo is the free electron mass. The spin factor is 
caused by the phase shift in the oscillation for up and 
down spins. When the Zeeman splitting becomes a half 
of the Landau level spacing, the amplitude of the oscilla- 
tion of the fundamental frequency (ja = 1) becomes zero. 
Eq. (1) is called the LK formulaB The experiments 
are performed in the multi energy-band system with the 
fixed electron number (iV) , where the chemical potential 
varies as a function of the magnetic field. However, the 
oscillation of the chemical potential is very small when 
the Fermi surface has a three-dimensional shape. Then 
we can apply eq. (1) or the superposition of eq. (1) with 
some frequencies even for the system with fixed electron 
number. 

On the other hand, if the system is two-dimensional or 
the interlayer coupling is much smaller than the Landau 
level spacings (this condition is satisfied in the quasi- 
two-dimensional system in the strong magnetic field), the 
oscillation of the chemical potential plays important roll. 
The magnetization for the fixed electron number is given 
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in the two-dimensional singlejhand system. The spin re- 
duction factor in this case iscllf 
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where [gm/2mo] is an integer part of gmjlrriQ. The spin 
reduction factor, Rp , is the periodic function of gm/2mo 
but is not the cosine function. 

If the magnetic field is tilted by 9 from the z-axis, the 
spin reduction factor is given by replacing g with gj cos 9 
in eqs. (2) and (4). The reduction factor for the pth 
harmonics, Rp 1 , is zero when cos9 = (pgm / mo) / (2n + 
1) with integer n. When cos6* = (gm/mo)/(2n + 1), 
Rp for odd p is zero. The tilting angle, 9, when Rp 1 
Rp becomes zero is called spin splitting zero condition 
which can be measured by tilting the magnetic field. 

When two energy bands exist in the two-dimensional 
system, the magnetization is not described by the su- 
perposition of the magnetization for thaauglo band-ieq. 
(3)) if the electroiijuimber is fixed BlPi^Hy>BE| 
One of the authorspElLP Harrison-et al.,0) Sandu et al.p 
Han et al.li^ and Fortin et alJlff treat the magnetic 
breakdown model as shown in Fig. 1(a), where the (3 
orbit is .obtained by the tunneling through the zone gap. 
NakancO and Alexandrov and BratkovskyEJ) have cal- 
culated the dHvA oscillation by using the independent 
two-band model, in which there exist two closed orbits 
(a and (3 oiiits^^ jshpaoL pi Fig. 1 (b). From these 
calculations]^' DMMbBB it has been shown that 
there exist frequencies in the magnetization (M(N, H)), 
which are forbidden in the scmiclassical theory. One of 
these forbidden frequencies corresponds to the area of the 
(3-ot orbit. On the other hand, the forbidden .osciHaij 



do not appear in M(/x, H) when (i is fixedl 
These forbidden oscillations are caused by the chemical 
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potential oscillation. 

The existence of the (3-a oscillation in the mag- 
netization has been observed experimentally in quasi- 
two-dipensional materials with two , nr. three energy 
bands^ k-(BEDIjTTF) 2 Cu(NCS^,BE3) a-(BEDT- 
TTF) 2 KHg(SCN) 4 EP and Sr 2 RuOz 




Fig. 1. (a) Fermi surface in the magnetic breakdown system, (b) 
Fermi surface in the two-band system. 



as -neglected to simplify the calculation 
o||ii|,|i2|,|iJ) j n orc Jer to study how the spin af- 
fects the "forbidden" [3-a oscillation and quantum mag- 
netic oscillations (a, (3, /3+a, 2a oscillations, etc.) in the 
two-dimensional magnetic breakdown systems, we have 
studied M(N—H) in the tight-binding model quantum- 
mechanically.EP This model can be applied to the quasi- 
two-dimensional organic conductors. We found that the 
Fourier transform amplitudes (FTAs) of the /3+a oscilla- 
tion is enhanced due to the Zeeman effect, which is quite 
different from the spin reduction factors, eq. (2) or eq. 
(4). However, as far as we know, the effect of the Zeeman 
term in the two-dimensional two-band model has never 
been discussed systematically. 

In this letter, we calculate M (N, H) at T = for the 
two-dimensional two-band model, where the magnetic 
breakdown is neglected. In this case neither chemical 
potential nor electron number for each band and each 
spin are fixed. Our purpose is to make clear how the os- 
cillation in M(N, H) are modulated by an electron spin. 



We study the model, whose energy at H = is given 



by 
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where e\ is the bottom energy of the energy band and 
i = a and (3 are the band index. The density of states 
for spin pi and the cyclotron effective mass m; are 



Pi = Cm,i,mi = ^Jmfm v 



(G) 



where C = S/2-kTi 2 and S is the real space area of the 
two-dimensional system. 

When the magnetic field (H Q ) is applied with the tilt- 
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Fig. 2. M(N, H) as a function of 1/H when g = 0, 1.0 and 2.6. 



ing angle 0, the z-component of the magnetic field is 
H = Hq cos 9. The energy spectrum is quantized as 

ei(H,n,a) =hu i (n+-+a{^)g i )+e b i , (7) 
2 mo 

where n is Landau index, uii = eH/miC, a = ±i and 
9i = cos 9. We take g a = gp and gi = g, although 
g a and gp are not necessarily the same due to the spin- 
orbit coupling. We set e — Ti = c = 1 and C = 1. 

When the electron number (JV) is fixed, we have to 
calculate the Helmholtz free energy F(N,H), which is 
given by 



F(N, H) = J2 PiKui a( H > n, <x), 



(8) 



where the summations are performed for the filled and 
partially filled energy levels. The electron number, AT, is 
decided by the chemical potential, p(0), at H = 0. The 
magnetization is obtained by 



M(N,H) = —dF(N, H)/dH. 



(9) 



We calculate M(N, H) by changing the value of g. When 
g = 0, M(N, H) are reduced to the spinless model.LJ 

We show M(N,H) for various values of g in Fig. 2, 
where m a /mo — 0.5,mp/mo — 1.0, e„ = 680,6^ = 
and ju(0) = 1000. We take C = 1 and 0.25 < 1/H < 
0.275, which corresponds to about 8 tesla and /x(0) ~ 0.2 
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eV. In these parameters, the area of the a and orbits 
corresponds to f a = 160 and fp — 1000. 

We calculate the FTAs of M(N, H), where there exit 
many combination frequencies, /3 ± a, (3 ± 2a, 2/3 ± a, 
2/3 ± 2a, in addition to a, /3, 2a and 2/3, which are shown 
in Fig. 3. 
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Fig. 3. The Fourier transform amplitudes of M(N, H) in Fig. 2. 



We show the FTAs of M(N, H) as a function of g in 
Figs. 4 and 5. We also show FTAs for other parameters 
(m a /mo = 0.65, mp/m$ = 1.0, e h a — 754, = and 
/Lt(0) = 1000) in Figs. 6 and 7. 




Fig. 4. The FTAs of some peaks (/„, f 2a ^ f@, hp, fp-a 
and f/3+ a ) in M(N,H) as a function of g for m a /mo = 
0.5, /mo = 1-0, e„ = 680 and = 0. 



It is clear from Fig. 4 that the ^-dependences of the 
FTAs are not described by the spin reduction factor for 
the single-band systems, eq. (4). For example, the FTAs 




Fig. 5. The FTAs of some peaks (f/3- 2a , f0+2a,j2/3- a , h/3+, 
hp— 2a and f20+2a) in M(N, H) as a function of g for m /rao 
0.5, 771,3 /m = 1.0, = 680 and A = 0. 



of the second harmonics oscillations (2a and 2/3 oscilla- 
tions) are not constant as a function of <?, which is not 
expected from eq. (4) for even p. 

The FTAs of the fundamental oscillations (a and /3 
oscillations) become zero when g = \{ra^/mi)(2n + 1), 
which can be seen in Figs. 4 and 6. This is the same as 
the condition of the spin splitting zero given by eq. (4) 
for p = 1 . 

We can see from Figs 4 and 6 that the /3+a and /3- 
a oscillations are suppressed when the FTAs of a or /J 
oscillations becomes zero. This is in contrast with the 
magnetic breakdown model (Fig. 4 in ref.B^). For exam- 
ple, in the magnetic breakdown model the FTA of the 
/3+a oscillation is large at g = 0.4, when the (3 oscillation 
is suppressed. The large amplitude of the /3+a oscilla- 
tion in the magnetic breakdown model might be under- 
stood as follows. The suppression of the /3 oscillation is 
caused by the cancellation due to the ir phase difference 
between the oscillations for the up spin and the down 
spin in the semiclassical picture. This dose not result in 
the cancellation of the /3+a oscillation in the magnetic 
breakdown systems because the /3+a oscillation corre- 
sponds to the larger orbit caused by the tunneling and 
the phase difference between the up and down spins in 
the /3+a orbit may not be n. In the two-band systems 
without the magnetic breakdown, the /3+a oscillation is 
not caused by the larger orbit but the chemical potential 
oscillation. As a result, the /3+a oscillation is suppressed 
in the two-band systems when a or (3 oscillation is sup- 
pressed. 

The FTA of the 2/3 + 2a oscillation is constant as a 
function of g, which can be seen in Figs. 5 and 7. 

By tilting the magnetic field the suppression of the 
peak of the /3 + a oscillations and the constant peak of 
the 2/3+2a oscillation as a function of g may be observed 
in Sr2RuC>4. __ . 

The Yamaji effectlij) should be also considered when 
the magnetic field is tilted if the Fermi surface has weak 
three-dimensionality. NakancclP shows that the /3 + a 
oscillations become large due to the Yamaji effect when a 
and /3 oscillations are enhanced in the spinless model. In 
quasi two-dimensional materials, the interplay between 
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the spin effect studied in this paper and the Yamaji effect 
should be considered. 

Ohmichi et alE3 have measured the FTAs of each os- 
cillations in the magnetoresistance in S^RuCU, where 
there is no suppression of the peak of the (3 ± a oscilla- 
tions due to the spin although the enhancement of each 
oscillations due to Yamaji effect is seen. Our theory for 



the magnetization cannot 
netoresistance experiment 
turn interference oscillatioi 



e compared with their mag- 
§_because the Stark quan- 



mm shou 



be considered in 
' have measured 



the magnetoresistance. Yoshida et al 
the angle-dependence of the FTA of the a oscillation. 
We expect that the effect of the spin and quasi-two- 
dimensionality can be observed in the angle-dependence 
of the FTAs of (3 ± a and higher harmonic oscillations of 
the dHvA oscillation in Sr 2 Ru04. 




Fig. 6. The FTAs of some peaks (f a , f 2a ^ f/3, f2/3, f/3- 
and f/3+a) in M(N,H) as a function of g for m a /mo 
0.65, mp /m = 1.0, e* = 754 and = 0. 




Fig. 7. The FTAs of some peaks (fg- 2a , f/3+2a,j20-a, hp+a, 
f20-2a and /2/3+2a) m M(N, H) as a function of g for m a /mo = 
0.65, mp/m = 1.0, e b a = 754 and e& = 0. 



In conclusion, we study the dHvA oscillation includ- 
ing the Zeeman effect for the two-dimensional two-band 
model, where the magnetic breakdown is neglected. We 
find the anomalous ^-dependences of the FTAs of (3+a, 
(3-a and 2f3+2a oscillations in M(N,H), where g = 
(grrii) / (2mo cos 9) . The (3 ± a oscillations are suppressed 



when the a or (3 oscillation disappears, and the FTA of 
the 2/3+2a oscillation is constant as a function of g. The 
effect of the spin on M(N,H) cannot be described by 
the spin reduction factor for the single-band (eq. (4)), 
whereas that in M(/x, H) is given by eq.(2). Even in the 
two-band model, however, the spin splitting zero condi- 
tion for fundamental frequencies (a and (3 oscillations) 
in M(N,H) is given by eq. (4). We expect that the in- 
dependences of the amplitudes of these oscillations may 
be observed in the experiment of tilting magnetic field in 
two-dimensional multi-band system such as Si'2Ru04. 
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